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Abstract 

We revisit the problem of interplay between the strong and the Coulomb interac- 
ts - 

tion in the charged-to-neutral yield ratio for BB and DD pairs near their respective 

thresholds in e + e~ annihilation. We consider here a realistic situation with a resonant 
interaction in the isospin 1 = channel and a nonresonant strong scattering amplitude 
in the 1 = 1 state. We find that the yield ratio has a smooth behavior depending on 
the scattering phase in the 1 = 1 channel. The same approach is also applicable to the 
KK production at the <^(1020) resonance, where the Coulomb effect in the charged-to- 
neutral yield ratio is generally sensitive to the scattering phases in both the isoscalar 
and the isovector channels. Furthermore, we apply the same approach to the treatment 
of the effect of the isotopic mass difference between the charged and neutral mesons 
and argue that the strong-scattering effects generally result in a modification to the 
pure kinematical effect of this mass difference. 
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1 Introduction 



The J PC = 1 resonances near the new flavor thresholds: T(AS), -0(3770), and 0(1020) are 
the well known sources in e + e~ experiments of pairs of the new-flavor mesons: respectively 
BB, DD, and KK. A number of experimental approaches depends on the knowledge of the 
relative yield of pairs of charged and neutral mesons: 

a(e + e" -> P+p-) 



pc/r. 



a(e+e- -> P°P°) 



where P stands for the pseudoscalar meson, i.e. B, D, or K, and dedicated measurements 
of such ratio have been done at the T(AS) resonance [1] at -0(3770) [2] and at 0(1020) [3]. 

The values of the ratio R c ^ n at all three discussed resonances are close to one due to these 
resonances being isotopic scalars, and it is the deviation of the discussed ratio from one that 
presents phenomenological interest. This deviation is generally contributed by the following 
factors: the isospin violation due to the Coulomb interaction between the charged mesons 
and due to the isotopic mass difference between charged and neutral mesons, and, in the 
case of the KK production at the 0(1020) resonance, a non-negligible nonresonant isovector 
production amplitude. The latter effect can be studied and described as the "tail of the p 
resonance", while the isospin breaking due to the mass difference is usually accounted for 
as a kinematical effect in the P wave production cross section factor p 3 , where p is the the 
cm. momentum of each of the mesons. The Coulomb effect has attracted a considerable 
theoretical attention. The expression for this effect in the ratio R c / n in the limit, where the 
resonance and the charged mesons are considered as point-like particles pf] has the simple 
textbook form: 

SR c/n = ^, (2) 

with a being the QED constant and v the velocity of each of the (charged) mesons in the 
cm. frame. However for the production of the real-life mesons the analysis is complicated 
by the charge form factors of the mesons[5], by the form factor in the vertex of interaction of 
the resonance with the meson pair [3, E] and generally by the strong interaction between the 
mesons [7J [8j [9] . In particular, it has been argued [8j [9] that the modification of the Coulomb 
effect by the strong (resonant) interaction between the mesons is quite significant. The 
previously considered picture of the strong interaction was however somewhat unrealistic. 
Namely, it has been assumed [8j [9] that the wave function in the 1=1 state of the meson 
pair is vanishing at short but finite distances, which would correspond to a singular behavior 
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of the strong interaction at finite distances. In this paper we derive the formulas for the 
Coulomb effect in the ratio R c l n under the standard assumption about the strong scattering 
amplitude in the channels with 1 = and 1=1. We find that in the case of the T(4S) and 
0(3770) resonances, where the heavy meson pairs are produced by the isotopically singlet 
electromagnetic current of the corresponding heavy quark, the strong-interaction effect in 
the Coulomb correction depends on the scattering phase 8\ in the 1 = 1 channel and is a 
smooth function of the energy across the resonance, while in the case of the Kaon production 
at and near the 0(1020) there is also a smooth dependence on the nonresonant part of the 
strong scattering phase 8q in the isoscalar channel inasmuch as there is a contribution of the 
isovector production amplitude at these energies. In either case we find that the behavior 
of the Coulomb effect is smooth on the scale of the resonance width, unlike the behavior 
previously found [8], [9] under less realistic assumptions. 

We further notice that essentially the same calculation can be applied to considering the 
effect on the ratio R c l n of the isotopic mass difference Am between the charged and neutral 
mesons, at least in the first order in Am, by considering the mass difference as a perturbation 
by a (constant) potential. In this way we find that the result coincides with the linear in Am 
term in the ratio of the kinematical factors p 3 only in the limit of vanishing strong scattering 
phase. Once the latter phase is taken into account, there arises a correction whose relative 
contribution is determined by the parameter (pa) with a being the characteristic range of 
the strong interaction. We therefore conclude that the conventionally used p 3 approximation 
for this effect may be somewhat applicable to the KK production at the 0(1020) resonance, 
where p ~ 120 MeV, but becomes quite questionable for the DD production at the -0(3770), 
where p ps 280 MeV. 

The strong-scattering phase in the P-wave state of mesons produced in e + e~ annihilation 
near the threshold is proportional to p 3 . We therefore expect the discussed effects of the 
strong interaction in the ratio R c ^ n to exhibit a measurable variation with energy. A mea- 
surement of this variation can thus provide an information on the strong scattering phases, 
which is not readily available by other means. 

The material in the paper is organized as follows. In Sec. 2 we consider the production of 
meson pairs by an isosinglet source and derive the formula for the correction to R c ^ n due to 
a generic isospin-violating interaction potential V(r) viewed as a perturbation. In Sec. 3 we 
generalize this treatment to the situation where the source is a coherent mixture of / = and 
1 = 1. The specific expressions corresponding to the Coulomb interaction and the isotopic 
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mass difference are considered in Sec. 4. Sec. 5 contains phenomenological estimates of the 
constraints on the parameters of the strong interaction between heavy mesons based on the 
currently available data [TJ [2] for BB and DD production. Finally, in Sec. 6 we summarize 
our results. 

2 General formulas for an isoscalar source 

We start with considering the behavior of the scattering wave functions of a meson-antimeson 
pair in the limit of exact isotopic symmetry, i.e. neglecting any Coulomb effects and the 
isotopic mass difference. We adopt the standard picture (see e.g. in the textbook [10]). 
where the strong interaction is confined within the range of distances r < a, so that beyond 
that range, at r > a the motion of the mesons is free. The two relevant independent solutions 
to the Schrodinger equation at r > a for the radial wave function in the P wave are the free 
outgoing wave 



and its complex-conjugate, f*(pr), describing the incoming wave. A general wave function 
of a pair of neutral mesons, (j) n {r) as well as of a pair of charged mesons, 4> c {f), in this region 
is a linear superposition of these two solutions. 

In the region of strong interaction, i.e. at r < a, the isotopic symmetry selects as 
independent channels the states with definite isospin, 1 = and 1=1, corresponding to 
the wave functions 0o — 0c + <f) n and = <p c — n . The detailed behavior of the 1 = and 
1 = 1 wave functions inside the strong interaction region is not important for the present 
treatment, and the important point is that the non-singular at r = 'inner' wave functions 
match at r = a particular linear superpositions of the incoming and outgoing waves (which 
superpositions in fact correspond to standing waves): 



where 6q and 5i are the strong scattering phases in respectively the isoscalar and isovector 
states. 

Consider now the production of meson pairs by a source localized inside the region of 
strong interaction, i.e. r < a, such as e.g. the electromagnetic current. The wave function 




(3) 




(4) 
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of the produced meson pairs at r < a is then determined by both the source and the 
strong interaction, and the relevant solution to the Schrodinger equation is chosen by the 
requirement that asymptotically at large distances, r — > oo, only an outgoing wave is present. 
Let us first consider the simple case where the relevant electromagnetic current is a pure 
isotopic singlet, which is the case for DD and BB pair production. Then in the limit of 
exact isotopic symmetry the outgoing waves for the 'n' and the 'c' channels have exactly the 
same amplitude, which for our present purpose can be chosen as one: 

/A(°)/y) = f(jpr) and (r) = f(pr) at r — > oo , (5) 



c 



where the superscript (0) stands for the approximation of exact isotopic symmetry. It can be 
noted that the approximation of the free motion beyond the region of the strong interaction 
in fact makes the expressions in Eq.Q applicable at all r > a, i.e. all the way down to the 
matching point r = a. It is helpful to notice for a later discussion that at the matching point 
the 1=1 wave function is vanishing while the 1 = function <f)^ contains only the outgoing 
wave. When continued into the strong interaction region, i.e. at r < a, the function O 
evolves into the solution determined by the strong interaction and the source. 

The isospin-violating effects of the Coulomb interaction and of the mass difference Am 
between the charged and neutral mesons can be generally described as being due to a presence 
of an extra potential V(r) in the 'c' channel beyond the region of the strong interaction: 
V = —a/r for the Coulomb interaction effect and a constant potential V = 2 Am describing 
the mass difference. In other words the wave function <p n of the 'n' channel is still determined 
at r > a by the radial Schrodinger equation for free P-wave motion_|, while the equation for 
the 'c' channel function <b r reads as 



dr 



2 +p 2 -mV(r)-^\ <f) c (r) = . (6) 



It is assumed throughout the present consideration that the isospin-breaking potential 
exists only at distances beyond the range of the strong interaction, i.e. that V(r) has support 
only at r > a. The justification for such treatment is that in the region of the strong force 



2 Clearly, in the considered here first order in the isospin violation only the difference of the interaction 
between the two channels is important, thus any such difference can be relegated to one channel, while 
keeping the other one unperturbed. Also, any effect of the mass difference in the kinetic term p 2 /m is of 
order v 2 /c 2 as compared to the discussed here effect of Am in the overal energy difference between the two 
channels, and is totally neglected in our treatment. 
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small isospin-violating effects are compared to the energy of the strong interaction, so that 
the contribution of any such effects arising at r < a is very small, while in the region r > a 
the relative contribution of the potential V(r) is determined by its ratio to the kinetic energy 
of the mesons, which is small near the threshold. 

It should be emphasized that although the interaction at distances r > a is present only 
in the 'c' channel, the wave functions in both channels are modified in comparison with those 
in Eq.fjSJ), as a result of the coupling between channels imposed by the boundary conditions 
at r = a. According to the setting of the problem of production of the meson pairs by a 
localized source, the appropriate modified functions are those containing at r — > oo only the 
outgoing waves 

c -(l + a;)/(pr), <j> n - (1 + y) f(pr) , (7) 

where the (complex) coefficients x and y arise due to the potential V, and are proportional to 
V in the considered here first order of perturbation theory. These coefficients determine the 
ratio of the production amplitudes: A c /A n = 1 + x — y, and the discussed here modification 
of the yield ratio: 

R c / n = l + 2Rex-2Rey . (8) 

The modified wave function in both channels is subject to two conditions: 
i: The channel with neutral mesons has only an outgoing wave at all r > a. In other words, 
the expression for 4> n (r) in Eq.(j7| is valid at all r down to r = a; 

ii: The wave function of the channel with isospin I = 1 at r < a should be proportional to 
the standing- wave solution matching the function xi i n Eq.fll]), since there is no source for 
the 1 = 1 state of the meson pairs. 

These two conditions are sufficient to fully determine the modified functions at r > a and 
thus to find the coefficients x and y. 

The first order in V(r) perturbation of the wave function in the channel with charged 
mesons is found in the standard way, using the P wave Green's function G + (r, r') satisfying 
the equation 

(J^+P 2 -^) G ^,r') = 8{r-r>), (9) 

and the condition that G+(r, r') contains only an outgoing wave when either of its arguments 
goes to infinity. The Green's function is constructed from two solutions of the homogeneous 
equation, i.e. from the functions f(pr) and f*(pr), as 

G+(r, r') = ^- lf(pr) f*(pr') 6{r - r') + f(pr') f*{pr) 6{r' - r)] , (10) 
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where 9 is the standard unit step function. The perturbation 5(j) c is then found as 

S(j) c (r) = m / G + (r,r') V(r') f{pr') dr' . (11) 



One readily finds from this explicit form of the solution that 5(fi c contains only the outgoing 
wave at asymptotic distances r — > oo: 



S(fi c I™ = -— f(pr) I V(r') \f(pr')\ 2 dr' , (12) 
so that the coefficient x is purely imaginary: 

X= ~Tvla V ( r ')\f(P r ')\ 2dr ' ( 13 ) 

and gives no contribution to the ratio of the production rates i? C//n described by Eq.fJHjjf. 

Consider now the matching of the wave functions at r = a. In this region of r one has 
r < r' in the integral in Eq. ljTT]) so that the correction in the 'c' channel has only an incoming 

wave: 

Wc(r) \ r ^ a = V r(pr) (14) 

with 

V = I V ( r ') [fiPr'f dr ' ■ ( 15 ) 
The wave functions <fio = <fi c + (fin and (fii = (fi c — (fin corresponding to the states with isospin 
1 = and 1 = 1 are then found as 

(fio \ r ^ a = 2/(pr) + y f(pr) +rj f*(pr) and fa \ r ^ a = rjf*(pr) - yf(pr) . (16) 

One can now apply the condition ii to determine the coefficient y. Indeed, the condition for 
the wave function (fii at r — > a to be proportional to f*(pr) + e 2lSl f(pr) requires y to be 
given by 

y = - V e 2 ^. (17) 

Upon substitution in Eq.® this yields 

2 



R c/n = 1 + - Im 
v 



oo 



e 2i5i / e 2i P r i + J_ ] V{r)dr 

J a V pT 



3 It can be noticed that the integral in Eq. (|13[) is divergent, which corresponds to the infrared-divergent 
behavior of the perturbation for the phase of the wave function, logarithmic for the Coulomb interaction 
and linear for a constant potential. This slight technical difficulty can be readily resolved, for our present 
purposes, by introducing an infrared regularizing factor exp(— Ar) in the potential and setting A — > in the 
end result. 
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3 Mixed isoscalar and isovector source 



The formula (|T8|) gives the general expression for the isospin-breaking effect in the considered 
yield ratio for the case where the mesons are produced by an isoscalar source. The presented 
consideration can also be extended to a situation where the source is a general coherent mix- 
ture of an isoscalar and isovector. The specific isotopic composition of the source determines 
the ratio of the coefficients of the amplitudes of the running outgoing waves in the 1 = 1 
and 1 = channels at the matching point r = a, which ratio we denote as A\/Aq, thus 
defining Ai and A as the production amplitudes in the respective channels (in the limit of 
exact isotopic symmetry). In this situation the generalization of the expressions in Eq.Q for 
radial wave functions in the 'outer' region r > a in the zeroth order in the isospin violation 
can be written as 

#(r) = (A + A,) f(pr) and 0<°>(r) = ( Aq - A,) f(pr) . (19) 

The isospin violation in the asymptotic form of these wave functions at r — > oo can then be 
parametrized, similarly to Eq.(jZD, by complex coefficients x and y as 

C - (A + A,) (1 + x) f(pr), n - (A Q - A,) (1 + y) f(pr) , (20) 

so that the yield ratio is found from 

j- 4, 2 

[l + 2Rex-2Rey). (21) 



R c/n A> + A x 2 



A -At 

The coefficient x, similarly to the previous discussion and the equation f[T51) . is purely 
imaginary and in fact does not contribute in Eq. (l2TI) . while the coefficient y is found from the 
appropriately modified conditions on the wave functions. Namely, the previously discussed 
condition i remains applicable, so that the asymptotic expression in Eq. fl2TJ|) for the 'n' 
channel function remains valid in the entire 'outer' region r > a down to the matching point 
r = a. In order to allow for the isovector component of the source the condition ii has to be 
modified as will be described few lines below. 

The perturbation by the potential V(r) of the 'c' channel wave function at the matching 
point r = a is readily found, similarly to Eq. tflll) . as 



with f] given by Eq. (|T5|) . 



Ti(A + A 1 )f*(pr) (22) 
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One can now write the expressions for the resulting 'outer' wave functions in the isotopic 
channels at the matching point: 



and 



<fo( r ) 



2 A f(pr) + v(A + At) f*(pr) + y (A - A 1 )f(pr) 



2A + y(A - A x ) - 7] {Ao + A x ) e 2i5 ° f(pr) + rj {A + A,) e*° X o(r 



JSo 



(23) 



2 Ax f(pr) +r](A + Ax) f*(pr) -y{A - A x )f(pr) 



2A 1 -y(A - A x ) - V (A + Ax) e 2 ^ f(pr) + v (A + Ax) Xl (r) 



i5i 



(24) 



with Xo an d Xi being the standing wave functions from Eq.(J3]) in the corresponding isotopic 
channels, which when evolved in the region of strong interaction contain no singularity at 
r = 0. The remaining parts in the latter expressions for the functions 0o an d 4>x describe 
the proper running outgoing waves. These parts, when continued down in r into the strong 
interaction region evolve to match the source at r < a. The ratio of the amplitudes of 
the isovector and the isoscalar running waves is determined by the isotopic composition of 
the source, and by the isotopically symmetric propagation through the strong-interaction 
region. Thus the ratio of the amplitudes of these waves at r = a does not depend on the 
isosp in-breaking effects at r > a and should be equal to Ax/A . Applying this condition to 
the isotopic wave functions given by the expressions fT2"3l and (jUj) , one finds the equation 
for the coefficient y: 

2 Ax - y (A - Ax) - V (A + Ax) e 2i _^ = Ax 

" A ' 



(25) 



2 Ao + y (Ao - Ax) ~v(A + Ax) 

This equation in fact replaces in this more general situation the previously discussed condi- 
tion ii, which condition and the ensuing result in Eq.( fl71) are readily recovered in the limit 
Ax/Aq = from Eq.M. 

Considering that both y and r] are of the first order in the potential V, it is sufficient to 
use the linear expansion of the equation ff25l) in y and rj, finding in this way the solution for 
y in the form 

Aq e 2lSl - Ax e 2i5 ° 



V = -"- A - Al 
and thus arriving at the final formula for the relative yield: 



R c/n 



Aq + Ax 



Aq-Ax 



1 + -Im 

v 



A e 2iSl - Axe 
Ao-Ax 



2i8 



7 2ipr 



1 + — V(r) dr 



(26) 



(27) 
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Given that Aq = \Aq \ e t5 ° and A\ = \Ax \ e %Sl , the amplitude-dependent factor in this formula 
can also be written in terms of the real ratio p = \Ai/Aq\ as 

A e 2iSl - Ai e 2i5 ° 2iS 1 - p e i( - 5a - 5 ^ 



A -A 1 



P*_ 

1 — pe 



-i(So—Si) 



(28) 



4 The Coulomb and the mass- difference effects 

The general formulas in Eq. (fl8|) and (|27|) can now be applied to a discussion of the specific 
isospin-breaking effects in the e + e~ production of meson pairs at and near the threshold 
resonances. We start with considering the effect of the Coulomb interaction. In a detailed 
treatment of this correction one should include the realistic form factors of the mesons, 
which cut off at short distances the difference in the electromagnetic interactions between 
the charged and neutral mesons. In the present discussion we replace for simplicity the 
gradual cutoff of the Coulomb interaction by an abrupt cutoff at an effective range r = a c , 
where generally a c > a0. The master integral with the Coulomb potential V(r) = —a/r in 
the equations fflBl and (J27j) then takes the form 



2ipr 



1 + — V(r) dr 



Ci(2pa, 



c sin2pa c 



cos 2pa c ( sin 2pa c 
pa c 

2(pa c f H2mc) 
where the integral sine and cosine are defined in the standard way: 





Si(2pa c 



(29) 



Si(z) 



dt 

smt — 

o t 



and Ci(z) 



00 dt 
cost — 
t 



and 7£ = 0.577 ... is the Euler's constant. The latter line in Eq. fl2"9l shows few first terms of 
the expansion of the integral in the parameter (pa c ). This expansion illustrates the behavior 
of the correction toward the threshold. For the purpose of this illustration one can consider 
first the simpler expression in Eq. ffl8|) . The imaginary part, which determines the discussed 
Coulomb effect in R c l n in the limit where there is no strong scattering, b\ — * 0, is not singular 



4 As previously mentioned, any extension of the isospin-breaking potential inside the strong interaction 
region can result only in very small corrections. 
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at pa c — > 0, and the textbook formula ([2]) is recovered in this limit. The real part of the 
integral in Eq. (129|) is singular at small pa c , but it multiplies in Eq. (jl8j) the factor sino~i. The 
P-wave scattering phase in its turn is proportional at small momenta to p 3 : 5i ~ {pa.) 3 , so 
that the overall contribution of the real part of the integral is not singular at the threshold 
either. Considering a more general expression for the Coulomb effect for the case of an 
isotopically mixed source, following from the equation ff27j) . one can readily arrive at the 
same conclusion that the singular in (pa c ) real part of the integral ( |29i) does not lead to 
an actual singularity, since it only enters the ratio P C//n multiplied by a combination of the 
phases 5q and 5i (cf. Eq. (1281) ) . each vanishing as p 3 toward the threshold. 

As previously mentioned, the effect of the isotopic mass difference corresponds to that of 
a constant potential V = 2 Am extending from the range of the strong interaction r = a to 
infinity. The master integral with such potential has the form 



2 sin 2 pa 

— cos 2 pa 




Am f 2 cos 2 pa 

+ sin 2 pa + i 



p { pa 
Am { 2 



p [pa 



pa 

2pa + 3i + {pa) 2 \ . (30) 



In the limit of vanishing strong scattering phases the mass correction to R c / n is determined 
by only the imaginary part of the integral, which in the limit of small pa thus yields 

vp E ' v ; 

where E is the total kinetic energy of the meson pair, and the found expression coincides 
with the linear in Am term in the expansion of the usually assumed ratio of the kinematical 
factors (p+/po) 3 - Clearly, in the more realistic case of presence of the strong scattering the 
real part of the integral in Eq. fl30|) also contributes and the simple kinematical approximation 
is generally invalidated. 



5 Phenomenological estimates 

In this section we discuss application of our formulas to interpreting the data on the charged 
to neutral meson yield ratio R c ^ n at the near-threshold resonances T(AS), ?/>(3770) and 
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0(1020). The purpose of this discussion is to illustrate the effect of the strong scatering on 
the isospin breaking corrections, and we use here the simplified picture of a abrupt cutoff 
of the Coulomb interaction and of the isotopic mass difference effects. Such simplification 
generally can be used as long as the parameter (pa) is not large. A detailed analysis should 
likely involve a model of a gradual cutoff, since the details of the transition become important 
at lager momenta. 

5.1 T(45) 

The simplest case for the study of the isospin breaking corrections in the relative production 
of heavy mesons is offered by the BB pair production near and at the T(45) resonance. 
Indeed, this process only is due to the purely isosinglet electromagnetic current of the b 
quarks, and the isotopic mass difference between the B mesons is very small: Am# = 
—0.33 ± 0.28 MeV [TT], so that any deviation of the ratio R c l n from one is essentially entirely 
due to the Coulomb interaction. On the other hand, the parameter a/v for the Coulomb 
effect in this case is the largest due to small velocity of the B mesons: at the energy of the 
T(4S) peak vb/c ~ 0.06. In particular, the numerical value in the expression §2§ is 0.19. 
The experimental data [1] however indicate a significantly smaller deviation of R c / n from one. 
The BaBar data with the smallest errors give R c ^ n = 1.006 ±0.036 ±0.031. Such behavior is 
likely a result of a combined effect of the meson and production vertex form factors [5], [6] and 
of the discussed here modification of the Coulomb correction by the strong scattering phase. 
These effects can in principle be separated and studied quantitatively by measuring the 
energy dependence of the ratio R c ^ n near the T(45) resonance. With the presently available 
data we can only use a simplified parametrization of the form factor effects by introducing 
an abrupt cutoff for the Coulomb interaction at r = a c > a and thereby estimate the likely 
regions in the (a c , Si) plane. Such estimate from the equations (TT8"]) and ( 1291) is shown in Fig.l 
as a one-sigma area, corresponding to the BaBar data with the statistical and systematic 
errors added in quadrature: R c/n = 1.006 ± .048. Clearly, more precise data from dedicated 
measurements of the ratio i? C//n are needed for a better understanding of the parameters of 
strong interaction between the B mesons. 
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0.2 0.4 0.6 0.8 1 1.2 

a c fm 



Figure 1: The one sigma area (shaded) in the (a c , 5\) plane corresponding to the BaBar data 
on the B + B~ /B°B° yield ratio at the T(4S) resonance. 

5.2 ^(3770) 

The largest isospin-breaking effect in the DD production at the ^(3770) is that due to the 
mass difference between the charged and the neutral D mesons: Amp = 4.78±0.10 MeV [llj . 
The most precise measurements of this process have been done [2] at the energy a/s = 
3773 MeV. At this energy the momentum of each charged D meson is p + = 254 MeV and 
that for a neutral D meson is po = 287 MeV. Thus the ratio of the kinematical factors 
(p+/po) 3 ~ 0.69 is significantly less than one. The Coulomb effect is somewhat smaller. 
Indeed, the velocity of a charged meson at this energy is v + /c = 0.135 and the expression 
(T2]) gives numerically 0.085. One can notice that if the kinematical and the Coulomb factors 
are combined in a straightforward way to estimate R c/n = {p+/pof [1 + 7ra/(2v+)] « 0.75, 
this would be in a very good agreement with the experimental number [2]: R°^ n = 0.776 ± 
0.024taoo6- Thus it is quite likely that at this particular energy there is a considerable 
cancelation between the strong-interaction effects in the yield ratio, and such cancelation by 
itself imposes constraints on the parameters of strong interaction between the D mesons, 
which constraints is interesting to analyze. 
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An analysis of the strong-interaction effects along the lines discussed in the present paper 
generally runs into two difficulties. One is that our approach is accurate only in the linear 
in Am approximation, while the actual effect of the isotopic mass difference between the 
D mesons is not very small. However, numerically, the first term in the expansion of the 
kinematical factor (Eq. (l3ip ) gives 0.67, which is quite close to the mentioned above value 
0.69, and it looks like the linear term gives a reasonable approximation. The other point 
is that the cutoff parameter a c for the Coulomb interaction at short distances does not 
necessarily coincide with the range parameter a used for the short-distance cutoff of the 
effect of the mass difference. However, as previously noted, the Coulomb effect is somewhat 
small at the energy of the ^(3770) resonance, and for the purpose of preliminary estimates 
we set a c = a in our numerical analysis. In order to allow for possible errors introduced by 
our approximations in comparing with the data, we linearly add a theoretical uncertainty 
of 0.03 units to the combined in quadrature statistical and experimental errors. Proceeding 
in this way we find that the only region in the (a, 5\) plane at a < 2 fm consistent with the 
CLEO-c data at one sigma level is the one shown in Fig. 2. 



4 0° 
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20° 
10° 


-10° 

0.2 0.4 0.6 0.8 1 

a fm 

Figure 2: The area (shaded) in the (a, 5i) plane corresponding to the CLEO-c data on the 
D + D~/D°D° yield ratio at the ^(3770) resonance. The uncertainty shown includes a one 
sigma experimental error with our estimate of the theoretical uncertainty added linearly. 
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It is interesting to compare the plots in the Figures 1 and 2. In the heavy quark limit 
applied to both b and c quarks the strong interaction between the heavy mesons should be 
the same, corresponding to the same range parameters a and a c . The scattering phase Si 
for these two systems is generally different due to different masses. However, provided there 
are no isovector 'molecular' bound states, the sign of the phase should be the same, with 
the absolute value of the phase for heavier B mesons being larger than for the D mesons. 
The comparison with the data for the D mesons favors small values of the range parameter, 
as indicated by Fig. 2. If one also assumes that a c ~ a for the B mesons, the short range 
of a c , according to Fig.l, is compatible with the B mesons dcitci £it a negative scattering 
phase 5i, which sign of 5\ is also in agreement with the D meson data. A negative sign of 5\ 
corresponds to a repulsion, which for the 1=1 state of heavy meson pairs can be expected 
on general grounds [12] . 

5.3 0(1020) 

We believe that the production of KK pairs in e + e~ annihilation at and near the 0(1020) 
resonance merits a separate analysis along the lines discussed in the present paper and using 
detailed data similar to those in Ref.[3]. As is known, this production receives a small but 
measurable nonresonant contribution from the isovector part of the electromagnetic current 
of the u and d quarks, which corresponds to an isotopically mixed source. Furthermore, it 
has been pointed out [13] that a detailed theoretical analysis of the K + K~ /K°K° yield ratio 
at the 0(1020) resonance produces a result which possibly is at a meaningful variance with 
the data. 

At present we limit ourselves to noticing that the formula in Eq. fT27|) . applicable in this 
situation, describes a smooth behavior of the considered isospin breaking effects across the 
resonance in the 1 = channel. Indeed, the 1 = scattering phase at energy E near the 
resonance energy E is given by the Breit-Wigner formula 

e ** = ^llgatfo (32) 
A + Z7 

where A = E — Eq, 5q is the nonresonant scattering phase in the isoscalar channel, and 7 is 
the width parameter. Both 5q and 7 are smooth functions of the energy proportional to p 3 at 
small momentum, and 7(-Eo) determines the resonance width T as 7 = T/2. The ratio of the 
isovector and isoscalar production amplitudes can then be parametrized near the resonance 
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as 



M _ A + ci(Si _~ 5o) 

Aq fl 



(33) 



where /i is a parameter with dimension of energy: /i ~ — m p . The amplitude ratio 
entering the correction factor in Eq. fT27|) can then be written in the form 



which manifestly shows that this ratio is a pure phase factor of a complex quantity slowly 
varying across the 0(1020) resonance. 

6 Summary 

We have considered the effects of the isospin breaking by the Coulomb interaction and by 
the isotopic mass difference in the relative yield i? C//n of pairs of charged and neutral mesons 
near threshold by a compact source, such as in the production of heavy mesons in e + e~ 
annihilation. These effects are modified by the strong interaction scattering phases. The 
general formula for a situation where the source is an arbitrary coherent mixture of an 
isoscalar and isovector is given by Eq. (f27I) . In particular, for a purely isoscalar source, which 
is the case for the e + e~ annihilation into DD and BB pairs the strong-interaction effect is 
determined by the scattering phase 5\ in the 1 = 1 channel (Eq.( fl8l) ). As a practical matter 
we find that under the standard assumptions about the strong scattering amplitudes in the 
near-threshold resonance region the ratio R c l n has a smooth behavior with energy showing 
no abnormal rapid variation on the scale of the resonance width. The energy dependence 
of this ratio is rather determined by the non- resonant scattering scattering phase(s). In 
the P-wave the phase 5\ is proportional to p 3 , so that a measurement of the behavior ratio 
R c ' n with energy can provide information on this phase, which is not readily accessible by 
other means. The behavior of the ratio R c ^ n at larger energies away from the threshold also 
depends on the details of the onset of the strong interaction between the heavy mesons at 
short distances and on the behavior of their electromagnetic form factors, and a study of 
this behavior can provide an insight into these properties of the heavy-light hadrons. 



A e 2iSl - At e 2i5 ° 

Aq — Ay 



//-(A -i 7) e~ i{Sl - 5o) 
/i- (A + i^) e +i ( Sl -So) 



(34) 
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